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We extend Birkhoff 's theorem for almost LRS-II vacuum spacetimes to show that the rigidity of 
spherical vacuum solutions of Einstein's field equations continues even in the perturbed scenario. 
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I. INTRODUCTION 

The core content of BirkhofF's theorem (l| is that any 
spherically symmetric solution of the vacuum field equa- 
tions has an extra symmetry: it must be either locally 
static or spatially homogeneous. This underlies the cru- 
cial importance in astrophysics of the Schwarzschild so- 
lution, as it means that the exterior metric of any exactly 
spherical star must be given by the Schwarzschild metric 
(it cannot be Minkowski spacetime if the mass is non- 
zero); and this also underlies the uniqueness results for 
non-rotating black holes. 

However it is an exact theorem that is only valid for ex- 
act spherical symmetry; but no real star is exactly spher- 
ically symmetric. So a key question is whether the result 
is approximately true for approximately spherically sym- 
metric vacuum solutions. We prove an "almost Birkhoff 
theorem" in this paper that shows this is indeed the 
case, so those results carry over to astrophysically re- 
alistic situations (such as the Solar System). There are 
of course many papers discussing perturbations of the 
Schwarzschild solution, but none appear to focus on this 
specific issue. It is in a sense an analogue of another 
important result, the "almost EGS theorem" [2] that 
generalizes from exact isotropy of radiation to approx- 
imate isotropy the crucial Ehlers-Geren-Sachs theorem 
Q, proving that isotropic cosmic background radiation 
everywhere in an expanding universe domain U implies 
a Robertson- Walker geometry in that domain. In both 
cases an exact mathematical result, depending on exact 
symmetry, is generalized to a more physically realistic 
result, depending on approximate symmetry. 

The almost-Birkhoff result will not be true if spacetime 
is not a vacuum (empty) spacetime, for the degrees of 
freedom available through a matter source generically in- 
validate the result, as is shown for example by the family 
of Lemaitre-Tolman-Bondi (LTB) models j^]. However it 
remains true for the electrovac case Q and solutions with 
a cosmological constant. It will also not be true for grav- 
itational theories with a scalar degree of freedom, such as 
scalar-tensor theories 0. The rigidity embodied in this 
property of the Einstein Field Equations is specific to 



vacuum General Relativity solutions, or solutions with 
a trace-free matter tensor. One should also note that 
(like the EGS result) the result is local: both Birkhoff's 
theorem, and our generalization of it, are independent of 
boundary conditions at infinity: they hold in local neigh- 
borhoods U . There are of course crucial differences be- 
tween the exactly spherical case and that of approximate 
symmetry: specifically, an almost spherically symmetric 
pulsating star can emit gravitational waves, but this is 
not possible in the exactly spherical case. However we 
show that in such a case if the intensity of the radiation 
is such as to leave the solution almost spherically sym- 
metric, the space time will remain almost static, which 
is defined in a way we will make precise later. 

We prove the result by using the 1+1+2 covariant per- 
turbation formalism 0,131 1 which developed from the 1+3 
covariant perturbation formalism Q. This enable us to 
prove the approximate result as a straightforward gener- 
alization of the exact result, which we prove first, using 
the 1+1+2 formalism. Actually we prove a small gener- 
alization of the standard Birkhoff result: it holds for all 
Class II Locally Rotational Symmetric (LRS) spacetimes 
[lOj (which include Schwarzschild as a special case). 



II. 1+1+2 COVARIANT SPLITTING OF 
SPACETUVIE 

The 1+3 covariant formalism 9] has proven to be a 
very useful technique in many aspects of relativistic cos- 
mology. In this approach first we define a timelike con- 
gruence by a timelike unit vector (u°'Ua = — 1). Then 
the spacetime is split in the form R^V where R denotes 
the timeline along and V is the tangent 3-space per- 
pendicular to Then any vector X"" can be projected 
on the 3-space by the projection tensor h'^ = Qfi + u°-ui,. 
The vector is used to define the covariant time deriva- 
tive (denoted by a dot) for any tensor T°' '^c..d along the 
observers' worldlines defined by 



rpa..h 



a..b 



(1) 



and the tensor hah is used to define the fully orthogonally 
projected covariant derivative D for any tensor T°'' '^c..d , 
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with total projection on all the free indices. 



(2) 
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In the (1+1+2) approach we further spHt the 3-space 
V, by introducing the spacehke unit vector e" orthogonal 
to so that 



eau" = , 
Then the projection tensor 



eae-" = 1. 



ga + Uav!' 



(3) 



2, (4) 



The hat-derivative is the derivative along the e° vector- 
field in the surfaces orthogonal to u". The 5 -derivative 
is the projected derivative onto the orthogonal 2-sheet, 
with the projection on every free index. We can now 
decompose the covariant derivative of e° in the direction 
orthogonal to into it's irreducible parts giving 



Daeb = CaOh + -(pNab + ^Eab + Cab , 



projects vectors onto the tangent 2-surfaces orthogonal 
to e*^ and u", which, following Q, we will refer to as 
'sheets\ Hence it is obvious that e"'Nab = = u°'Nab- In 
(1+3) approach any second rank symmetric 4-tensor can 
be split into a scalar along a 3-vector and a projected 
symmetric trace free (PSTF) 3-tensor. In (1+1+2) slic- 
ing, we can take this split further by splitting the 3-vector 
and PSTF 3-tensor with respect to e°. Any 3-vector, tj:"' , 
can be irreducibly split into a component along e° and a 
sheet component 5'", orthogonal to e'' i.e. 

^a = ^e" + , * = V'^ea , = iVVb ■ (5) 

A similar decomposition can be done for PSTF 3-tensor, 
■006, which can be split into scalar (along e"), 2- vector 
and 2-tensor part as follows: 



where 



Cab 



e Dc^a — e-a , 
s: a 

^ ^ab X 

■^S OaSb , 

^{a^b} ■ 



(12) 

(13) 
(14) 

(15) 
(16) 



We see that along the spatial direction Ca, (j> represents 
the expansion of the sheets Cab is the shear of (i.e. the 
distortion of the sheet) and a" its acceleration. We can 
also interpret ^ as the vorticity associated with e" so that 
it is a representation of the "twisting" or rotation of the 
sheet. The other derivative of e° is its change along u". 



Ca -^Wq 



(17) 



V'ab = V'{ab) 

where 



1 



*a 
*ab 



e"eVab = 
iVa'e^Vbc 



N^'^ab 



2^'(aeb) + *ab , (6) where we have A = e'^Ua and aa = NacC'^- Also we can 
write the (1+3) kinematical variables and Weyl tensor as 
follows 



= N\aNbf - -^NabN"' i^cd , (7) 



e = 



u 



hybu'' 

Ae'' + A'' , 



(18) 
(19) 
(20) 



and the curly brackets denote the PSTF part of a tensor 
with respect to e". 
We also have 



1 



(Jab = S ( CaCb - -iVab ) + '^^ia^b) + ^ab , (21) 



h{ab} = , A^(ab) — -e(aeb) — Nab 



2, 

^hab 



(8) 



Eat 
Hab 



S l^eaCb - 
T-L [ CaCb 



The sheet carries a natural 2-volume element, the alter- 
nating Levi-Civita 2-tensor: 



]^Nai}j + 2£(,eb) + £ab , (22) 
^Nab^ + 2H(,eb) + Hab ■ (23) 



.d 



Vdabce U 



(9) 



where Sabc is the 3-space permutation symbol the vol- 
ume element of the 3-space and "qabcd is the space-time 
permutator or the 4-volume. With these definitions it 
follows that any 1+3 quantity can be locally split in the 
1+1+2 setting into only three types of objects: scalars, 
2-vectors in the sheet, and PSTF 2-tensors (also defined 
on the sheet). 

Now apart from the ' time^ (dot) derivative of an object 
(scalar, vector or tensor) which is the derivative along 
the timelike congruence u", we now introduce two new 
derivatives, which defines, for any object i}a...b'^"''^'- 



where Eab and Hab are the electric and magnetic part of 
the Weyl tensor respectively. Therefore the key variables 
of the 1+1+2 formalism are 

'^ab,£ab,T~iab,C,ab] ■ (24) 

These variables (scalars , 2-vectors and PSTF 2-tensors) 
form an irreducible set and completely describe a vacuum 
spacetime. In terms of these variables the full covariant 
derivatives of e" and m° are 



Sfipa. 



c.d 



c.d 



e^Df'iJja..b 



C.d 



(10) 



Naf...Nb'Nh'.-N,^NfW,4,f„g'-^ .(11) 



VaGfc = -AUaUb - UaCtb + ( S + -6 ) CaUb 

+ (Sa - eac^") Ub + CaOb + ^(j)Nab 
+^£ab + Cab , 



(25) 
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VaUfa = -Ua {Acb + Ab) + CaBb ( ^9 + S 



+ VLSab + Safe . (26) 



For the complete set of evolution equations, propaga- 
tion equations, mixed equations and constraints for the 
above irreducible set of variables please see equations (48- 
81) of [1]. Also we have the following commutation re- 
lations for the different derivatives of any scalar ?/', 

= -A'4> + {^Q + ^)i:+{Y.a+eac^''-aa)5''i^ , (27) 

SaSbi' - SbSai' = ieabi^i' " C^) + 2a[Jb]i^ (28) 

From the above two relations it is clear that the 2-sheet 
is a genuine two surface (rather than just a collection 
of tangent planes) , in the sense that the commutator of 
the time and hat derivative do not depend on any sheet 
component and also the sheet derivatives commute, if 
and only if + Sac^'^ — Q.a ~ and 51 = ^ = a" = 0. 



III. BIRKHOFF THEOREM FOR VACUUM 
LRS-II SPACETIMES 

Locally Rotationally Symmetric (LRS) spacetimes 
posses a continuous isotropy group at each point and 
hence a multi-transitive isometry group acting on the 
spacetime manifold These spacetimes exhibit lo- 

cally (at each point) a unique preferred spatial direction, 
covariantly defined. Since LRS spacetimes are defined to 
be isotropic about a preferred direction, this allows for 
the vanishing of all orthogonal l-|-l-f 2 vectors and ten- 
sors, such that there are no preferred directions in the 
sheet. Then, all the non-zero 1-1-1-1-2 variables are covari- 
antly defined scalars. A subclass of the LRS spacetimes, 
called LRS-II, contains all the LRS spacetimes that are 
rotation free. As consequence in LRS-II spacetimes the 
variables fi, ^ and H are identically zero and the vari- 
ables {A, 8, (j), fully characterize the kinematics of 
the vacuum spacetime. The propagation and evolution 
equations of these variables are: 



£ = - 



1 9 

2^^ 



-e-s 



A- 



s-^e =- ^0 + 2( ie-^E) -£, 



£ ^ 



(29) 

(30) 
(31) 

(32) 

(33) 
(34) 



A-e =- {A 



\A + -e'^ 



3v.2 

2^ 



(35) 



Since the vorticity vanishes, the unit vector field 
is hypersurface-orthogonal to the spacelike 3-surfaces 
whose intrinsic curvature can be calculated from the 
Gauss equation for u" that is generally given as (ll| : 



^^^Rabcd = iRabcd)± — KacKbd + KbcKad 



(36) 



where ^^^Rabcd is the 3-curvature tensor^ _L means pro- 
jection with hab on all indices and Kab is the extrinsic 
curvature. Also we note that in case of LRS-II space- 
times the sheets at each point mesh together to form 

2- surfaces. The Gauss equation for together with the 

3- Ricci identities determine the 3-Ricci curvature tensor 
of the spacelike 3-surfaces orthogonal to u°' to be 



Rab 









eaSb- 



K 



This gives the 3-Ricci-scalar as 



'R 



3 

4^ 



K 



Nab . (37) 



(38) 



where K is the Gaussian curvature of the sheet, '^Rab = 
KNab ■ From this equation and (|29p an expression for K 
is obtained in the form [llj 



K = -£ 



(39) 



From (|29ll34p . the evolution and propagation equations 
of K can be determined as 



K 



-9 - S K 



K 



(40) 
(41) 



From equations (|3Tj) , (|34l) , (gOl) and ([IT]) we get an inter- 
esting geometrical result for vacuum LRS-II spacetime 



-(t)K. 



£ = CK^/"^. 



(42) 



That is, the H-1+2 scalar of the electric part of the Weyl 
tensor is always proportional to a power of the Gaussian 
curvature of the 2-sheet. The proportionality constant C 
sets up a scale in the problem. We can immediately see 
that for Minkowski spacetime C = 0. 

To covariantly investigate the geometry of the vacuum 
LRS-II spacetime, let us try to solve the Killing equation 
for a Killing vector of the form — ^''"a + ^Sa, where 
\l/ and $ are scalars. The Killing equation gives 

Va(«'U6-f-$eb)+Vb(^'Ua + $ea)=0. (43) 

Using equations (|25l) and (PB)) . and multiplying the 
Killing equation by u'^e'', e°'e^ and iV"'' we get the 
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following differential equations and constraints: 



solution 



w - 


- = 


L) , 


(44) 


^ _ $ _ + 




0;, 


(45) 


l> + *(-e 

3 


+ s) = 


0, 


(46) 




f $0 = 


. 


(47) 



Now we know ^aS." = + If C° is timelike (that is 
faf° < 0), then because of the arbitrariness in choosing 
the vector w", we can always make $ = 0. On the other 
hand, if is spacelike (that is ^a^" > 0), we can make 
^' = 0. 

Let us first assume that i^" is timelike and $ = 0. In 
that case we know that the solution of equations pi)) 
and (|45|) always exists while the constraints (|46)) and 
(|T7)l together imply that in general, (for a non trivial 
^'), 8 = S = 0. Thus the expansion and shear of an unit 
vector field along the timelike Killing vector vanishes. In 
this case the spacetime is static. Now if is spacelike 
and \1/ = 0, solution of equations ps)) and always 
exists and the constraints pij) and (|Tf| together imply 
that in general, (for a non trivial $), (j) ^ A = 0. 
From the LRS-II equations (|29p - ((35)) . we then imme- 
diately see that the spatial derivatives of all quantity 
vanish and hence the spacetime is spatially homogeneous. 

In other words, wc can say: There always exists a 
Killing vector in the local [u, e] plane for a vacuum 
LRS-II spacetime. If the Killing vector is timelike then 
the spacetime is locally static, and if the Killing vector is 
spacelike the spacetime is locally spatially homogeneous. 

In fact in the first case, when = E = 0, we have 
K — 0. Furthermore if choose coordinates to make the 
Gaussian curvature 'iiT' of the spherical sheets propor- 
tional to the inverse square of the radius co-ordinate 'r', 
(such that this coordinate becomes the area radius of the 
sheets), then this geometrically relates the 'hat' deriva- 
tive with the radial co-ordinate 'r'. As we have already 
seen, K = —4>K, where the hat derivative, defined in 
terms of the derivative with respect to the co-ordinate 
'r', depends on the specific choice of e" (orthogonal to 
M° and the sheet). If we choose the 'radial' co-ordinate 
as the area radius of the spherical sheets, then from pH)) 
and (|14p the hat derivative of any scalar M becomes 



M 



1 dM 



(48) 



for a static spacetime. If the spacetime is not static then 
there is also a dot derivative in the RHS of equation . 
Details are given in [TTj |. 

Now solving equations (P^ - ([55]) . we get the unique 



2ra 
r 

E = — 



A 



m 



2 m 
r 



(49) 
(50) 



Here the constant m is the constant of integration. Solv- 
ing for the metric components using the definition of 
these geometrical quantities we get [11| 



ds^ = -{l-'^\dt^ 



dr^ 



(1 



2m ^ 



^dn'^ 



(51) 



which is the metric of a static Schwarzschild exterior. 
In the second case, when <j) — A ~ 0, we can choose 
_ ^ 2m _ -j"^a .^^j^gj.g jg a coustaut and then solving 
(1^ - we get the unique solution 

3m — 2t 2 3m — t 



e = 



t^t{2m~t) 
_2m 



ity/t{2m-t)' 



-4 



Again solving for the metric components we get 
dt^ 



ds^ 



(¥-1) 



2m 



1 dr^ +t^dn'^. 



(52) 
(53) 

(54) 



which is a part of the Schwarzschild solution inside the 
Schwarzschild radius. 

Thus we have proved the (local) Birkhoff Theorem: 

Any solution of Einstein's equations in empty space, 
which is of the class LRS-II in an open set S , is locally 
equivalent to part of maximally extended Schwarzschild 
solution is S. Also it is interesting to note that the 
modulus of the proportionality constant in equation (j42p , 
which sets a scale in the problem, is exactly equal to the 
Schwarzschild radius. 



IV. ALMOST BIRKHOFF THEOREM 

As we have already seen, for a spherically symmetric 
spacetime, the 1-1-1-1-2 scalar of the electric part of 
the Weyl tensor is always proportional to the (3/2)*'' 
power of the Gaussian curvature of the 2-sheet, with the 
proportionality constant defining a scale in the problem. 
Let us now perturb this geometry and define the notion 
of an almost spherically symmetric spacetime in the 
following way: 

Any spacetime that admits a local 1+1+2 splitting 
at every point such that the magnitude of all 2-vectors 
on the sheet, the sheet gradient of scalars (defined by 
\/'0a'0°y'; magnitude of all PSTF 2-tensors on the 
sheet, and the sheet derivative of 2-vectors ( defined by 
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\/^pabi^) dt any given point are either zero or much 
smaller than the scale defined by the modulus of the 
proportionality constant in equation ^J^, is called an 
almost spherically symmetric spacetime. 

We would like to emphasize here that though 
Minkowski spacetime belongs to the set of LRS-II, in the 
above definition of the perturbed spacetime we exclude 
the Minkowski background, as in that case the scale is 
identically zero. As we have seen from equations ([?7|) 
and the sheet will be a genuine two surface if and 
only if the commutator of the time and hat derivative do 
not depend on any sheet component and also the sheet 
derivatives commute. In the perturbed scenario we will 
require the sheet to be an almost genuine 2-surface in the 
sense that the commutator of the time and hat derivative 
almost do not depend on any sheet component and the 
sheet derivatives almost commute; this will follow from 
our definition of almost spherically symmetric. In that 
case the scalars and ^ would be of the same order of 
smallness as the other vectors and PSTF 2-tensors on the 
sheet. Also using the constraint equation 

San^' + eabS"^'' = {2A~(t>)n-3^^+eabC"^c+'H , (55) 

we see the scalar Ti also is of the same order of smallness. 
Hence the set of 1+1+2 variables 

8"','H°',T,ab,Sab,'HabXab] ■ (56) 

are all of 0{e) with respect to the invariant scale. Using 
equations (48-81) of Q, we can get the propagation and 
evolution equations of these small quantities. We now 
list all the equations up to the first order (that is up to 
order 'e'). The time evolution equations of ^ and C{a6} 
are as follows: 

e = (is - ^e) ^+{A~^<j))n + ieab^'a^ + , (57) 

C{a6} = (iS-ie)Ca6+(-4-i0)S„b 

+S{aab} - £c{a^b} ■ (58) 

The Vorticity evolution equations: 

tl = iea&<5M'' + AS, + VI {Y. - le) , (59) 

na + ^SabA' = -(|0+iS)r!a 

+ieab [-Aa'> + S'^A - ^M''] -(60) 
Shear evolution equations: 

^{ab} = S^aAb} + ACab - ip + ^ab - Sab 161) 



Evolution equation for e^: 

aa-da = - {^(p + A) aa + {p + Y^) {Aa - aa) 

+ (i0-^) (S„+e„f,r!^) -e,^,H^ (63) 
Electric Weyl evolution: 

£a + I'ieabn — jEab^'^'H + ^EbcS^'H'^a ~ f^^a 

+ l£eab^'' - if «a + (|S ~~ B) £a 

-{\4> + A)eabn\ (64) 

£{ab} ~ £c{afi-b} — ~£c{aS'^T~(-b} ~ ^afc ^ (^ + f^) ^ab 

+ (i0 + 2A) e.^aU.f . (65) 
Magnetic Weyl evolution: 

n = -SabS'^E'' - 3^8 + (6* + |S) H , (66) 

^^a — ^EabS'^ — (f S — 9^ Ha ^ ^^EabA^ 

+ (i0 + A) EabE'' - leabS^'E , (67) 

Hiab} + Scia^k}" = +l^£c{aC,^ - + 2^) ec{a'^^6}'= 

- (^+|S]) Hah + Ecia-^^} • (68) 

In the above equation all the zeroth order quantities 
are background quantities. If the background is static 
with 6 = E = and the time derivative all the back- 
ground quantities are zero, we can easily see that the 
time derivatives of the first order quantities at a given 
point is of the same order of smallness as themselves. 
Hence the first order quantities still remains "small" as 
the time evolves. 

Similarly we can write the spatial propagation equa- 
tion of all the first order quantities up to 0{e). The 
propagation equations of f and C{a6} are: 

i = -<^e +{\e + Y)n + IsabS'^a'' , (69) 

C{ab} = -4'Cab+ 5{aab} + {\0 + Y)Yjab~ ^ab ■{^^) 

Shear divergence: 

+ {lc^ + 2A)eab^''-5''^ab, (71) 



ta-\Aa = \5aA+{A~\(f))Aa-{le+\Y.)Ya 

+ \Aaa - lY.aa - £a (62) 



S{ab} = (^{aSf,} - ec{a5''^b} " \<t>'^ab 

+ |SCa6 - ec{aHb} ■ C''^) 
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Vorticity divergence equation: 

Ci = -Sail'' + {A-4>)n . (73) 

Electric Weyl Divergence: 

£a = ^Sa£-Sb£ab-^£aa-^(l)£a. (74) 

Magnetic Weyl divergence: 

H = -SaW - - 3£n , (75) 

Ha = l25an~S''nab-l£eab^''+l£^a 

+ I^Sab£' - ■ (76) 

Hence we see that if the background is spatially homo- 
geneous with (j) = A = Q and the 'hat' derivative all the 
background quantities are zero, we can easily see that the 
'hat' derivatives of the first order quantities at a given 
point are of the same order of smallness as themselves. 
Hence the first order quantities still remain "small" along 
the spatial direction. In both the cases of a static back- 
ground and a spatially homogeneous background the re- 
sultant set of equations are the perturbed LRS-H equa- 
tions, (that is equations (|29l) - (|35l) with ©(e) terms added 
to each). 

Again trying to solve the Killing equation (|43p for a 
Killing vector of the form = '^Ua + ^Ca, using equa- 
tion ((26| . (|25|) and multiplying the Killing equation by 
^t"^t^ ^i'^e^ e°e^ iV'^^ 7V°^t^ N^e^ and iV°^d. we get 
the following differential equations and constraints: 





= 0, 


(77) 


^-$-*yt + $(E + ie) 
3 


= , 


(78) 


$ + $(ie + i]) 

o 


= 0, 


(79) 


2 

'i'(-e - s) + 

3 


= 0, 


(80) 




= 0, 


(81) 


(5c$ + $ac + Z^'Sc 


= 0, 


(82) 




= . 


(83) 



Now we see that for both timelike ($ ~ 0) or spacelike 
(^' — 0) vectors, all the above equations are not com- 
pletely solved in general unless the first order quantities 
appearing in the equations above are exactly equal 
to zero. This special case corresponds to static but 
distorted black holes in the presence of matter outside 
black hole, for example when an accretion disk occurs 
[l^ . If the distribution of the matter outside black 
hole is axisymmetric, then the vacuum metric outside 
the matter is described by Weyl Solution. However as 
we proved that these first order quantities generically 
remain 0{e) both in space and time, we can see that a 
timelike vector with (0 = E = 0) or a spacelike vector 
with {4> = A = Q) almost solves the Killing equations. 



Therefore we can say: 

For an almost spherically symmetric vacuum space- 
time there always exists a vector in the local [u, e] plane 
which almost solves the Killing equations. If this vector 
is timelike then the spacetime is locally almost static, 
and if the Killing vector is spacelike the spacetime is 
locally almost spatially homogeneous. 

Also as we have seen that in this case the resultant 
set of equations are the perturbed LRS-II equations, 
(that is equations (1^ - with 0{e) terms added to 
each), and the perturbations locally remain small both 
in space and time, a part of the maximally extended 
almost-Schwarzschild solution will then solve the field 
equations locally. 

Thus we have proved the (local) Almost BirkhofT 

Theorem: Any solution of Einstein's equations in 
empty space, which is almost spherically symmetric in 
an open set S , is locally almost equivalent to part of a 
maximally extended Schwarzschild solution in S. 

Note that we do not consider perturbations across the 
horizon: our result holds for any open set S that does 
not intersect the horizon in the background spacetime. 
The result almost certainly holds true across the horizon 
also, but that case needs separate consideration. 

The above result can be immediately generalized in the 
presence of a cosmological constant. In that case an 'al- 
most' spherically symmteric solution in an open set S, is 
locally almost equivalent to part of a maximally extended 
Schwarzschild deSitter/anti-deSitter solution in S. Also 
the result holds for an almost spherically symmetric elec- 
tric charge distribution with no spin or magnetic dipole. 
In this case we have to use the energy momentum tensor 
of the electromagnetic field in vacuum with the magnetic 
part being equal to zero. Proceeding exactly in a similar 
manner one can then show that the solution of the per- 
turbed field equations will be almost equivalent to a part 
of maximally extended Reissner-Nordstrom spacetime. 



V. DISCUSSION 

The rigidity of spherical vacuum solutions of the EFE, 
as enshrined in Birkhoff's theorem, is maintained in the 
perturbed case: almost spherical symmetry implies al- 
most static. This is an important reason for the stability 
of the solar system, and of black hole spacetimes. 

Though there are many discussions in the literatures 
on the stability of Schwarzschild solution in General Rela- 
tivity, for example [l3| ; most of them deal with a specific 
sector of the maximally extended Schwarzschild mani- 
fold, namely the static exterior part. In this paper we 
have established in a compact and completely different 
way, an aspect of the stability of the static sector of 
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the complete manifold: as long as the solution remains 
almost spherically symmetric, it remains almost static, 
with a similar result for the spatially homogeneous sec- 
tor. Furthermore our result, being local, does not depend 
on specific boundary conditions used for solving the per- 
turbation equations. Hence it does not depend on the 
global topology of the spacetime, and brings out covari- 
antly the rigidity and uniqueness of the almost-spherical 



vacuum solutions of Einstein's field equations. 

This rigidity has interesting implications for the issue 
of how a universe made up of locally spherically symmet- 
ric objects imbedded in vacuum regions is able to expand, 
given that BirkhofF's theorem tells us the local spacetime 
domains have to be static. A two-mass exact solution il- 
lustrating this situation is given in p4| ; the present paper 
suggests the results given there are stable. 
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